Decoherence of domains and defects at phase transitions 
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In this further letter on the onset of classical behaviour in field theory due to a phase transition, we 
show that it can be phrased easily in terms of the decoherence functional, without having to use the 
master equation. To demonstrate this, we consider the decohering effects due to the displacement 
of domain boundaries, with implications for the displacement of defects, in general. We see that 
decoherence arises so quickly in this event, that it is negligible in comparison to decoherence due to 
field fluctuations in the way defined in our previous papers. 
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The standard big bang cosmological model of the early 
universe, with its period of rapid cooling, gives a strong 
likelihood of phase transitions, with concomitant symme- 
try breaking. This paper is a further paper in a sequence 
by ourselves and collaborators P, in which we ex- 
plore the way in which such phase transitions naturally 
take us from a quantum to classical description of the 
universe. 

That (continuous) transitions can lead rapidly to clas- 
sical behaviour is not surprising. Classical behaviour has 
two attributes, i) Classical correlations: By this is meant 
that the Wigner function(al) W[tt, cj>] peaks on classical 
phase-space trajectories, ii) Diagonalisation of the deco- 
herence functional, whose role is to describe consistent 
histories. 

Continuous transitions supply both ingredients. 
Firstly, the field ordering after such a transition is due 
to the growth in amplitude of unstable long-wavelength 
modes, which arise automatically from unstable maxima 
in the potential. From the papers of Guth and Pi [|| 
onwards, it has been appreciated that unstable modes 
lead to correlations through squeezing, and we shall not 
consider it further. Secondly, we understand diagonali- 
sation to be an almost inevitable consequence of coarse- 
raining. The stable short-wavelength modes of the field 
together with the other fields with which it interacts, 
form an environment whose coarse-graining enforces di- 
agonalisation and makes the long-wavelength modes de- 
cohere. It is how this is implemented that is the basis of 
this letter. 

We stress that there is more than one way to formu- 
late diagonalisation. In our earlier papers 

II 1 1 we 

required the density matrix p(t) itself to become diago- 
nal, rather than the decoherence functional. Whichever 
approach we adopt, there is an issue as to which field 
basis we attempt to enforce diagonalisation, which can 
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only be approximately achieved. For an infinite degree 
of freedom system field, ideally we should think of diag- 
onalisation functionally. In practice this is impossible to 
achieve, and we are forced to adopt a piecemeal approach 
in which we make 'mini-superspace' approximations in 
which a finite number of degrees of freedom are isolated 
as the most significant. On comparing these, the relevant 
ones are those which decohere last. Since we are look- 
ing for the onset of classical behaviour, we take adjacent 
classical solutions with which to require no quantum in- 
terference. With these caveats we shall see later that the 
two approaches lead to the same results. In each case 
a probabilistic description is obtained, but the approach 
given here permits easier calculation. 

Since phase transitions take place in a finite time, 
causality guarantees that correlation lengths remain fi- 
nite. In our previous work, and this, we use the formation 
of domains after a transition to characterise the onset of 
classicality. To recapitulate, we showed [l|, [U, [H, by us- 
ing the master equation for the reduced density matrix, 
that the environment renders the long- wavelength modes 
of the order parameter field classical at early times, by 
or before the transition is complete. In particular, those 
modes on the scale of the domain size will have deco- 
hered, even though the modes on the scale of domain 
boundary thickness do not. In this paper we recreate that 
result rather more simply and show that a parallel result 
holds due to small displacements of domain boundaries, 
in a coarse-graining that is insensitive to their positions. 

This latter result has another consequence. If the sym- 
metry breaking permits nontrivial homotopy groups the 
frustration of the order parameter fields is resolved by 
the creation of topological defects to mediate between 
the different ground states 0, @. Since defects are, in 
principle, observable, they provide an excellent experi- 
mental tool for determining how phase transitions occur. 
For the simple theory that we shall consider here, that of 
a real scalar field with double-well potential, the domain 
boundaries are the defects (domain walls) and we can 
view our results as decoherence as induced by small dis- 
placements of defects. The generalization to a complex 
field cf> is straightforward, and has been considered else- 
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where Q. This gives more substance to our preliminary 
attempts to show that vortex defects are also classical by 
the time of their production Q. 

As in our earlier work, we restrict ourselves to flat 
space-time. The extension to non-trivial metrics is 
straightforward in principle 

We now consider the case of a real quantum field <p 
with double-well potential in detail. As we have said, 
the field ordering after the transition begins is due to the 
growth of long- wavelength modes. For these modes the 
environment consists of the short-wavelength modes of 
the field, together with all the other fields with which 
6 inevitably interacts in the absence of selection rules 
9, 1(1 . The inclusion of explicit environment fields is 
both a reflection of the fact that a scalar field in isolation 
is physically unrealistic, as well as providing us with a 
systematic approximation scheme To be specific, the 
simplest classical action with scalar and environmental 
fields Xa is 

s[(f>,x] = s syst {4>} + s env [x] + s qu [<p,x], (i) 

where (with fi 2 , m 2 > ) 

SsysM = J d A x «w + \v?4> 2 - ^> 4 } , 



and the most relevant interactions between system and 
environment are of the biquadratic form 

N 

s^x] = -E i J d ' x ( 2 ) 

Even if there were no external x fields with a quadratic 
interaction of kind of Eq.©, the interaction between 
short and long-wavelength modes of the </>-field can be 
recast, in part, in this form, showing that such a term is 
obligatory. 

Although the system field </> can never avoid the de- 
cohering environment of its own short- wavelength modes 
0, to demonstrate the effect of an environment we first 
consider the case in which the environment is taken to be 
composed only of the fields x&- Since environments have 
a cumulative effect on the onset of classical behaviour, 
the inclusion of a further component of the environment 
reduces the time it takes for the system to behave classi- 
cally. Thus it makes sense to include the environment one 
part after another, since we can derive an upper bound 
on that time at each step. 

We have shown elsewhere [l|, [H, H| that, in order to 
make our calculations as robust as possible, we need a sig- 
nificant part of the environment to have a strong impact 
upon the system-field, but not vice-versa. The simplest 
way to implement this is to take a large number N 3> 1 of 



scalar Xa fields with comparable masses m a ~ fi weakly 
coupled to the <p, with A, g & <C 1 (for details see Ref.[l[). 
Thus, at any step, there are N weakly coupled environ- 
mental fields influencing the system field, but only one 
weakly self-coupled system field to back-react upon the 
explicit environment. 

For one-loop consistency it is sufficient, at order of 
magnitude level, to take identical g a = g/yN. Further, 
at the same order of magnitude level, we take g ~ A. This 
is very different from the more usual large-N 0(N + 1)- 
invariant theory with one (/>-field and TV Xa fields, dom- 
inated by the 0(l/N) (x 2 ) 2 interactions, that has been 
the standard way to proceed for a closed system fllj ] . 
With our choice there are no direct x 4 interactions, and 
the indirect ones, mediated by <f> loops, are depressed by a 
factor g/yN . In this way the effect of the external envi- 
ronment qualitatively acts as a proxy for the effect of the 
internal environment provided by the short-wavelength 
modes of the c/>-field, but in a more calculable way. 

We shall assume that the initial states of the system 
and environment are both thermal, at a high tempera- 
ture To > T c . We then imagine a change in the global 
environment (e.g. expansion in the early universe) that 
can be characterised by a change in temperature from To 
to Tf < T c . That is, we do not attribute the transition 
to the effects of the environment-fields. 

Given our thermal initial conditions it is not the case 
that the full density matrix has <f> and x fields uncorre- 
cted initially, since it is the interactions between them 
that leads to the restoration of symmetry at high temper- 
atures. Rather, on incorporating the hard thermal loop 
'tadpole' diagrams of the x (and </>) fields in the 4> mass 
term leads to the effective action for </> quasiparticles, 

s*m = J d± x |ia„^ - \ m \{T^ 2 - 

(3) 

where, in a mean-field approximation, m^(To) = —fi 2 (l — 
T 2 /T c 2 ) for T ps T c . As a result, we can take an initial 
factorised density matrix at temperature To of the form 
p[To] = p^po] x/3 x [Xb], where p^po] is determined by the 
quadratic part of [0] and /5 x po] by Senv[Xa]- That 
is, the many Xa fields have a large effect on 0, but the 
0-field has negligible effect on the Xa- 

Provided the change in temperature is not too slow the 
exponential instabilities of the 0-ficld grow so fast that 
the field has populated the degenerate vacua well before 
the temperature has dropped significantly below T c [l2l |. 
Since the temperature T c has no particular significance 
for the environment fields, for these early times we can 
keep the temperature of the environment fixed at T x ps 
T c (our calculations are only at the level of orders of 
magnitude). Meanwhile, for simplicity the Xa masses are 
fixed at the common value m ~ p. 

It is sufficient for our purposes here to take an instan- 
taneous quench. Slower quenches make the analytic cal- 
culations very much more difficult, without changing the 
qualitative nature of the results [jj . 
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The notion of consistent histories provides an alterna- 
tive approach to classicality to trying to solve the master 
equation, as we have done previously. Quantum evolu- 
tion can be considered as a coherent superposition of fine- 
grained histories. Since we need to be able to distinguish 
different classical system-field configurations evolving af- 
ter the transition, we work in the field-configuration ba- 
sis. If one defines the c-numbcr field <p(x) as specifying a 
fine-grained history, the quantum amplitude for that his- 
tory is ^f[4>] ~ e lS ^ (we work in units in which H = 1). 

In the quantum open system approach that we have 
adopted here, we are concerned with coarse-grained his- 
tories 



[a] = / V<P e 



ism, 



(4) 



where a[<p] is the filter function that defines the coarse- 
graining. In the first instance this filtering corresponds 
to tracing over the Xa degrees of freedom. 

^From this we define the dccohcrcncc functional for 
two coarse-grained histories as 



V\a + ,a 



V4> + V(f>~ e 



"])«+ 



[4>+]a 



(5) 

T>[a + ,a } does not factorise because the histories (fx* 1 
are not independent; they must assume identical values 
on a spacclike surface in the far future. Decoherence 
means physically that the different coarse-graining his- 
tories making up the full quantum evolution acquire in- 
dividual reality, and may therefore be assigned definite 
probabilities in the classical sense. 

A necessary and sufficient condition for the validity 
of the sum rules of probability theory (i.e. no quantum 
interference terms) is (ill ] 



ReV 



0. 



(6) 



when a + ^ a (although in most cases the stronger 
holds (3). Such histories are 



condition T>\a + , of 
consistent (15j [ . 

For our particular application, we wish to consider as 
a single coarse-grained history all those fine-grained ones 
where the full field <fi remains close to a prescribed clas- 
sical field configuration <j) c \. The filter function takes the 
form 



a c \[4>] 



VJ e^^-^ctciLJ]. 



(7) 



In the general case, a[(j>] is a smooth function (we exclude 
the case a[(f>] = const, where there is no coarse-graining 
at all). Using 



J<j) = / d i xJ(x)4>(x), 



(8) 



V[a + ,a-} = j VJ+VJ- jw[J + ,J-)-{^ 
x a + [J + ]a-*[J-}, 

where 



(9) 



a iW[J+,J-] _ 



- p i(s[4> + ]-s[4>-]+J+<p+-J-<p-) 



_ (10) 

is the closed- path-time generating functional [lq . 

In principle, we can examine adjacent general classi- 
cal solutions for their consistency but, in practice, it is 
simplest to restrict ourselves to particular solutions </>^, 
according to the nature of the decoherence that we are 
studying. Initially, as we said earlier, we have made a 
de facto separation into the order parameter field </> and 
its explicit environmental fields Xa whereby, in a saddle- 
point approximation over J, 



F\ 



'cl> VcU 



(11) 



F[(f) + ,(f>~] is the Feynman- Vernon [171 | influence func- 
tional (IF) (see Ref. [J, Q for details). The influence 
functional is written in terms of the influence action 

A [<f>tv<t>ci\ as 



F\<tti,4G] = exp{iA[(/>+ fa}}. 



As a result, 



exp{-Im5A[</)+,0 cl ]}, 



(12) 



(13) 



where SA is the contribution to the action due to the 
environment. 

^From this viewpoint, once we have chosen the classical 
solutions of interest, adjacent histories become consistent 
at the time to, for which 



1 w Im SA\ t 



(14) 



As we are considering weak coupling with the envi- 
ronment fields, we may expand the influence functional 
F[(j> + , up to second non-trivial order in coupling 
strengths for large N. Higher terms are depressed by 
powers of N. The general form of the influence action is 
then d [3 

6A[^+ t 4T] = {(5 int [0 + ,x^])o-<5i„t[r,Xa]>o} 



-i{(S in t[(f> + , xt]Sint[<t> ,Xa])o 
(Sint , xt])o(Sint , Xa ])o} 



(15) 



,-{<S? nt [<T,Xa]>0 " [(^i„t[0-,X a -])o] 2 }. 



we may write the decoherence functional between two 
classical histories as 



As a further step we make a separation of the </>-ficld it- 
self into its long wavelength modes with |A;| < fx, which 
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determine the domains, and the short wavelength modes 
(|fc| > fi) which act as their implicit environment. How- 
ever, in calculating tu this latter step only serves to fur- 
ther reduce decoherence times which, as we shall see, are 
already short enough. 

In reality, even if we ignore this step for the purpose of 
bounding tjj there is not a unique bound, since it depends 
on the class of classical solutions considered. However, 
in practice all sensible choices seem, qualitatively, to give 
the same upper bound on try- The reasons for this are 
the following. Firstly, for the long wavelength system 
modes the spatial profile of the classical solution is not 
particularly relevant, since it is its exponential growth 
in time that sets the scale for the onset of classical be- 
haviour. Thus, not only is to insensitive to wavelength 
for small-k modes, but it only depends logarithmically 
on the parameters of the theory. Further, with the com- 
mon mass-scale fi that we assume here, the Compton 
wavelength fi , both determines the rate of exponen- 
tial growth and provides the natural distance scale over 
which we do not wish to discriminate between classical 
solutions. Different choices lead to differences in fitn of 
order unity, significantly smaller than fitjy itself, which 
we ignore. For weak couplings it is relatively easy to 
compute the upper bound on this time due to the inter- 
actions 5i n t[^, x) °f ©• We shall now see that, in our 
particular model, it is in general shorter than the time t* 
at which the transition is complete, defined in terms of 
the system field as the time for which 



V 



2 = 6fi 2 /\ 



(16) 



where the average is taken over the system field. That is, 
in (|16|i we ignore the contributions from the dominant en- 
vironment given by the short wavelength (|fc| > fi) modes 
of <f>. This renders (4> 2 )f finite. (More explicitly, in the 
language of our earlier papers based on the master equa- 
tion, e.g. [l|, (4> 2 )f = Tr{p r (f> 2 }, where p Y is the reduced 
density matrix obtained on tracing out the environment.) 
Any ambiguities are again of order unity in fitrj, the time- 
lag for non-linear effects to occur, in comparison to the 
exponential runaway of the free field [l9l |. 

Since the effect of the environment is to induce damp- 
ing, the classical behaviour of the field is expressed 
through the classical Langevin stochastic equations that 
it satisfies [2(1 HJ]. We are assuming that coupling is 
sufficiently strong for the system not to recohere after t* 

M- 

For the biquadratic coupling of Eq.([2]), the IF is given 
by 



Re5A = 9 — 



d 4 x 



d 4 y A(x)K(x-y)Z(y), 



ImSA 
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d 4 x / d 4 y A(x)N(x,y)A(y), 



where K(x — y) = lmG 2 i _ + (x,y)9(y° — x°) is the dissi- 
pation kernel and N(x — y) = KeG 2 , + (x, y) is the noise 



(diffusion) kernel. G++ is the relevant closcd-timc-path 
correlator of the x-field at temperature To [l[. We have 
defined A = \{4>+ 2 - 0~ 2 ) and £ = ±(0+ 2 + 4>~ 2 ). 
We look for classical solutions of the form 



(/)cl(x,s) = f(s,t)®(x), 

where, in principle, f(s, t) satisfies /(0, t) = <f>i and 
f(t,t) = 4>{ and gives the space-field configuration. 

We begin by showing that, using (TTI|) ,we recreate the 
results obtained previously in Refs.p], on using the 
more complicated master equation, in which the field is 
spread through space, and decoherence is due to differ- 
ent field amplitudes. In anticipation that to < t*, it 
is sufficient to restrict ourselves to the initial Gaussian 
free field evolution with negative (mass) 2 . In fact, for 
idealised Langevin equations, this can be a good approx- 
imation for domain formation into the non-linear regime 
[l2T ]. We look for classical classical fields that, after a 
sudden quench, have the form Q 

<p cl (s, x) = e^ s (/)f cos(fco x) cos(fco y) cos(k z), (17) 

where 4>f is the final field configuration. This is a single 
mode approximation to a regular chequer-board domain 
structure. Shorter wavelengths can be introduced with- 
out altering the result significantly. The reader is referred 
to [l[ for more details. For an instantaneous quench, 
we will use the late time behaviour of the longest wave- 
lengths (fco = 0), (f>^(s,x) ~ e^ s <j)f. The exponential 
factor, as always, arises from the growth of the unstable 
long- wavelength modes. 
Thus, lmSA[^ c 



^ci> ^ci] takes the form 



Im 5 A = 



^XlllA 2 
64 f 



dk 



(fc 2 + /x 2 ) 2 



1 + e^ 1 - e 2 ^ cos (2^/fc 2 + p?) 



(k 2 



(18) 



-2/i 2 ) 

and T r is the critical 



where A f = l/2(0+ 2 - cf)' 2 ), 
environmental-temperature. As we noted in previous 
publications, the volume V is due to the fact we are 
considering field configurations spread over all space. V 
is interpreted as the minimal volume inside which there 
are no coherent superpositions of macroscopically distin- 
guishable states for the field. Later, we shall consider 
localized configurations where this factor does not ap- 
pear. 

After assuming fit 1, the integral in momenta can 
be done analytically obtaining, 

I m ^^ 2 ^^»A t y. (19) 



256 



/'■' 



With this expression at hand, we are able to evaluate the 
decoherence time try for amplitude variation as 



fit D 



16/^5 



gT c A f F3 7ry / (3-2 v / 2) 



(20) 
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Using a conservative value for the volume, V 
0(pi~ 3 ), we get 



fj,t D 



2 



16^ 3 
gT c Ai 



(21) 



We re- write last expression in terms of Af = 0A/2, with 
<j) = <j>f + (f)^ and A = <fi^ — <j)~^ . At the completion of 
the transition cp 2 ~ r? 2 ~ A , and we will adopt, at time 
tD, 4> 2 ~ C(/i 2 a/A). A < a < 1 is to be determined self- 
consistently from the condition that, at time to, 
arj 2 . We have shown in 
determined as a ~ n/T c , 
we do not discriminate between field amplitudes which 
differ by 0(/j)), where characterizes the thickness 
of domain boundaries (walls) as the field settles into its 
ground-state values. Additionally, for simplicity we take 
the couplings g ~ A. Therefore, we obtain an upper 
bound onto, 



Ref. that the value of a is 
We also set A ~ 2/x (i.e. 



fj,t D 



2 \T cV ^ 



(22) 



which exactly coincides with the result found in Ref. 2] 
from the master equation for a sudden quench but, in the 
present case, using the decoherence functional approach. 

If we now trace over short-wavelength modes, as in 
but for unstable modes, we would get a further term in 
Im 8 A qualitatively similar to (|19[) . which will serve to 
preserve to < t* by making to even smaller. We shall not 
consider this implicit environment in subsequent analysis. 



For comparison, we find t* 

by 0,1 



for which (d> 2 ) 



lit* 



2 



rj , given 



(23) 



whereby /i 1 < to < t* , with 



fit* - \xto 



Hi 



> i, 



(24) 



for weak enough coupling, or high enough initial temper- 
atures. 

Having seen that the decoherence functional approach 
gives identical conclusions to the solution of master equa- 
tions, the new work in this letter is to evaluate the IF 
for different field configurations which, from the point of 
view of the Master equation, would be much more tax- 
ing analytically. Instead of the classical solutions used 
before 0, in which the field is spread in an infinite 
chequer-board through space, here we are concerned with 
a different field configuration; a localized domain wall. 

The orientation of such a wall is irrelevant, as was the 
orientation of the chequer-board solution used earlier. 
Most simply, we consider classical domain wall solutions 
(for the fco = mode) of the form 



(frfe 1 * 8 tanh(^x), 



(25) 



which link adjacent domains. 

Our main new result is to determine the decoher- 
ence induced by a small displacement in the domain 
wall, by evaluating the influence action for the classi- 
cal field configurations (/> ± (x,s) = 77/(5, t)$ ± (x), where 
$ ± (a;) = $(2; ± <5/2). We will consider 5 as a small dis- 
placement in the position of the wall, and consequently 
we expand the classical solution (or, equivalently, A(cc)) 
in powers of fi5, up to linear order 



A(s,x) a \i S r] 2 e 2Als tanh(fix) sech 2 (/j,x) . 

Doing this, the imaginary part of the influence action 
can be written as (after integrating over time and assum- 
ing fit > 1), 



2 T 2 4r2 2 

l mSA a 9 t c l \? e 4 ^ 



d A k / d A p 



64(2tt) 6 

e -i{p+k){x-y) 



d x / dry 



(fc 2 + /i 2 )(p 2 + M 2 ) 
tanh(/^a;) sech 2 (fj,x) tanh(/i?/) sech 2 (/ly) 



(1//C 2 + n 2 + yjf- + /i 2 ) 2 + A^ 2 



These integrations can be exactly evaluated (in part 
analytically and in part numerically), to give 



ImSA a 0.2 



(7 2 T 2 ?7 4 L 2 J 2 e 4 ^ 4 



1024 



(26) 



where L 2 is a surface term, analogous to the volume V in 
the chequer-board analysis. This coefficient comes from 
the fact that we are considering a one dimensional kink 
solution embedded in a three dimensional space. The 
L 2 coefficient comes from the "free surface" or wall in 
two directions of the three-dimensional kink. Considering 
very conservative values (such as to obtain upper bounds 
for to), we can set L = C(/j _1 ) as the minimum length 
scale in which there could be not coherent superpositions 
of macroscopic states of the field. 

If the decoherence time due to displacements is to, let 
us suppose that [ito > 1- Then to can be evaluated from 
the last equation, and its order of magnitude is given by 



/ito 



1 



In 



74/^ 



gT c r) 2 L5 



(27) 



Superficially, this result looks similar to previous results, 
where the difference between the field amplitudes has 
been replaced by the scale of displacement suffered by 
the domain boundary wall, 6. The main difference is in 
the power of 77 inside the log. If one takes 6 = JpL , wc 
may bound to as 



ut n ~ — In I — — 

' 2 \ 7 T C 



(28) 
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Since T c ^> fi (in fact one can show T 2 //^ 2 ~ 24/A 
this result puts a bound on the possible value of 7 
in order to have filo > 1, i.e. 7 < n/T c <C 1, and hence 
8 <C M • ^ 7 i s larger, at very early times in our model 
t ~ /Lt _1 , we get from Eq. (f26|) Im <5A > 1 immediately. 
Thus, the system behaves classically from this earliest 
permissable time, fito = C(l)- 

Since the size of the core of the domain wall is ©(/i" 1 ), 
we take 8 ~ 0(fi~ 1 ) (i.e. 7 = 0(1)) as the minimum 
length scale in which there could be not coherent super- 
positions of macroscopic states of the field. With deco- 
herence occurring at the earliest possible time it follows 
automatically that 

fit* - fiI D > 0. (29) 

Even if we demand decoherence at scales a fraction of the 
domain wall core thickness (i.e. 7 < 1) we can ensure 
that 

fit* - fit D ~ 1 In 1 1^ J > 0. (30) 

This result implies that decoherence due to the displace- 
ment of the boundary is a very early time event Id < tjj. 

The outcome of this analysis is that configurations of 
field domains with displaced boundaries are less impor- 
tant for decoherence than fixed domains in which there 
are field fluctuations, which decohere later. 

What is more appealing is, if one perform all the calcu- 
lation in 1+1 dimensions, the coefficient L appears only 
in our first example (i.e. the plane- wave); not for the 
kink-like solutions. 

In fact, rather than consider a whole chequer-board of 
domains as in [l], [2], H[ , we can restrict ourselves to two 
adjacent domains with boundary given by (|25p . Con- 
sider the simple case where the difference between walls 
is given by <fif = <fi^~ — e, with e a small fluctuation around 
the final value of the field configuration. Wc assume that 
the final field configuration is (f>^ = yfck], with a the 
self-consistent coefficient shown earlier and e = O(fi). 



Perhaps surprisingly, we recreate the result of ([24]) ex- 
actly. Yet again, we have decoherent behaviour before 
the transition is completed. 

In addition, the field possesses classical correlations at 
early times by virtue of the quasi-Gaussian nature of the 
regime @, [23| to give a fully classical picture by time t* . 

Finally, we return to the other view of domain bound- 
aries mentioned in the introduction, as topological de- 
fects. What does this analysis have to say about the 
classical behaviour of defects like vortices, whose separa- 
tion also measures the size of domains in simple circum- 
stances [24|? Treating vortices as having their cores as 
line zeroes (in the same way that the cores of our domain 
walls are sheet zeroes) we have shown elsewhere that the 
mechanism for the production of classical vortices has 
several parts 0]. Yet again, the environment renders 
the long-wavelength modes of the order parameter field 
classical at early times, by or before the transition is 
complete. In particular, those on the scale of the separa- 
tion of the line-zeroes that will characterise the classical 
domains will have decohered by the time the transition 
is complete, even though the field modes on the scale 
of classical vortex thickness do not decohere ever. For 
line-zeroes to mature into vortex cores the field needs 
to have an energy profile commensurate with the vortex 
solutions to the ordinary classical Euler-Lagrange equa- 
tions. This requires that the field fluctuations are peaked 
around long- wavelengths, to avoid fluctuations causing 
wiggles in the cores and creating small vortex loops, a 
related condition satisfied in our models. The resultant 
density of line-zeroes can already be inferred in the linear 
regime, whose topological charges are well-defined even 
though close inspection of their interior structure is not 
permitted classically. What the analysis of this letter 
shows is that we expect decoherence due to vortex dis- 
placement to be irrelevant in comparison to decoherence 
due to field fluctuations. 
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